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1. 


INTRODUCTION 


Plans  for  using  the  space  shuttle  orbiter  vehicle  to  launch  a 
number  of  payloads  have  been  under  recent  investigation.  It  is  proposed 
that  the  shuttle  be  launched  from  the  Eastern  Test  Range  inserting  the 
orbiter  into  a park  orbit. 

An  interim  upper  stage  (IUS)  which  is  to  be  deployed  from  the 
orbiter  and  used  to  achieve  the  desired  final  orbit  has  been  of  interest 
to  the  Air  Force  as  well  as  a number  of  private  contractors. 

The  park  orbit  was  defined  by  a number  of  considerations  to  be 
a circular  orbit  with  altitude  of  160  n mi  and  inclination  37.4°.  It  was 
desired  that  the  final  orbit  have  an  inclination  of  63.  4°,  an  argument 
of  perigee  of  270°,  and  a period  of  twelve  hours.  The  specific  twelve 
hour  orbit  of  interest  has  an  apogee  altitude  of  21450  n mi  and  a perigee 
altitude  of  350  n mi.  It  became  apparent  that  most  of  the  proposed  IUS 
configurations  could  not  perform  the  orbit  transfer  from  park  orbit  to 
final  orbit  using  two  burns,  and  a till  yield  some  desired  payload  values. 
Consequently  a study  was  initiated  to  investigate  the  feasibility  of  an 
optimal  three  burn  orbit  transfer. 

This  report  presents  the  results  of  the  optimal  three  burn  orbit 
transfer  investigation.  Initially  the  definition  of  a simplified  orbit 
transfer  simulation,  using  Keplerian  orbit  transfers  and  impulsive 
velocity  inc rements,  is  given.  After  discussing  the  use  of  a nonlinear 
programming  algorithm  to  solve  this  simplified  problem,  the  results 
for  a family  of  final  orbits  with  differ  jnt  inclinations  are  presented. 

The  accuracy  of  the  simplified  simulation  has  been  verified  by  comparison 
with  a detailed  trajectory  simulation. 
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THE  ORBIT  TRANSFER  PROBLEM 


a)  Three  Burn  Impulsive  Velocity  Orbit  Transfer 

Let  us  begin  with  the  definition  of  a simplified  orbit 
transfer  simulation. 

1)  Specify  the  initial  state:  the  position  vector  y(tQ),  the 
velocity  vector  y(tQ),  and  the  initial  time  tQ.  We  assume  that  the 
state  is  specified  using  an  earth-centered  inertial  (ECI)  coordinate 
system,  and  is  a point  on  the  160  n mi.  circular  park  orbit  with 
inclination  37.4°.  Furthermore  we  shall  assume  that  a spherical 
earth  model  describes  the  geopotential  function. 

2)  Coast  in  the  orbit  until  time  = tj  = t0  + 4tj,  where  4t^  is 
the  length  of  the  first  coast  in  seconds.  Essentially  the  coa-  t is 
simulated  by  integrating  the  equations  of  motion,  y = -jilyN  y, 
from  tQ  to  t^.  The  result  of  the  integration  is  the  state  vector 
y(tj) , y(ti)  and  because  of  the  spherical  geometry  assumption  this 
integration  can  be  performed  quite  rapidly.  The  complete  algorithm 
for  two-body  motion  in  space  used  in  the  simulation  is  described 

In  Escobal 

3)  Simulate  the  first  burn,  by  adding  to  the  inertial  velocity 
vector  y(t^)  the  velocity  increment  defined  by  (dVj,  Jflj,  Ail'i). 

The  pitch  increment  and  the  yaw  increment  A are  defined 

with  respect  to  the  lpcal  inertial  velocity  as  shown  in  Figure  1,  and 
are  applied  as  a yaw -pitch- roll  sequence. 


hV 


Figure  i.  Definition  of  Yaw  Increment,  and  Pitch 

Increment,  AS 
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4)  Coast  in  the  first  transfer  orbit  as  defined  by  the  new 

state  vector  until  time  = = i + At2. 

5)  Simulate  the  second  burn,  i.  e.,  add  to  the  inertial  velocity 
vector  the  increment  defined  by  (AV2,  A02,  A^2). 

6)  Coast  in  the  second  transfer  orbit  until  time  = *3  = t2  + A t2 

7)  Simulate  the  third  burn  - add  the  increment  defined  by 

(Av3,  A03, 

8)  End  of  Simulation  - Final  Orbit. 

It  is  important  to  note  that  the  coast  times  Att  and  the  velocity 
increments  (AV^,  A 6^,  A4>^)  i = 1,2,3,  in  conjunction  with  the  initial 
conditions  completely  determine  the  final  orbit.  The  entire  sequence  of 

2 

events  was  simulated  using  the  Generalized  Trajectory  Simulation  (GTS) 

b)  The  Optimization  Problem 

The  sequence  of  steps  required  to  compute  the  final  orbit 
have  been  outlined.  Explicitly  one  must  choose  the  twelve  variables 

x = <At1(  At2,  At3,  AVj.  AV2,  AV3,  A6lt  A^,  A62,  A^,  AQj.A^) 


such  that  the  objective  function 

f(x)  - AVj  + AV2  + AV3  (1) 

is  minimized  and  the  constraints 

c^x)  = 1^  - 21,450  = 0.  (2) 

c2(x)  = Up  - 3 50  = 0.  (3) 

c3(x)  = C0p  - 270  = 0.  (4) 

c4(x)  = i - 63.  4 * 0.  (5) 

c5(x)  = 75,000  - h^tj)  * 0 (6) 

c6(x)  = 7 5,000  - ha(t2)  * 0 (7) 


are  satisfied,  where  h is  the  apogee  altitude  (n  mi),  h the 

a 1 p 

perigee  altitude  (n  mi),  Wp  the  argument  of  perigee  (deg),  and  i the 
inclination  (deg). 


>■ 
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The  first  two  constraints  on  apogee  and  perigee  altitude  and 
are  representative  of  the  family  of  twelve  hour  period  orbits  required. 
The  third  constraint  requires  that  the  argument  of  perigee,  to  be  270° 
and  the  fourth  constraint  defines  the  final  orbit  inclination.  The  last 
two  constraints  establish  limits  on  the  apogee  altitudes  of  the  first  and 
second  transfer  orbits  respectively,  and  are  included  to  prevent  escape 
orbits  during  the  transfer.  The  7 5000  n mi  limit  is  somewhat  arbitrary 
and  its  significance  will  be  discussed  later. 


3.  THE  NONLINEAR  PROGRAMMING  PROBLEM 

The  previous  section  outlined  an  example  of  a nonlinear  programming 
problem.  Stated  concisfly  the  general  problem  is  to  determine  the  n-vector 
x = (x  | , x2...*xn)  that  minimizes  (maximizes)  the  objective  function 

f(x)  = f(Xj,  x2,  ....  xn)  (8) 

subject  to  the  equality  constraints 

c ^(x)  = 0 i = 1 . . . k (9) 

and  the  inequality  constraints 

c^(x)  * 0 i = (k  + 1),  . . . , m.  (10) 

Define  the  Lagrangian 

L(x,  A)  = f(x)  + c^(x)  A (11) 

where  A is  the  m-vector  of  Lagrange  multipliers.  The  Kuhn-Tucker 

jJj  $ 

necessary  conditions  require  that  at  the  optimum  point  (x  , A ) 

7L(x*.  A*)  = g(x*)  + G(x*)  A*  = 0 (12) 
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where  g(x)  is  the  gradient  vector  of  f(x),  and  the  n x m Jacobian  matrix 
is  defined  by 

f6cl  6cml 


G(x)  = [Vcj,  . . . , V ] 


v)t  ) $1 

Also  at  (x  , X ) 


n*  • 

n 


\ tj  (x  ) =0 


where  A.  • s 0 


i = 1,  ....  m 


i = (k+  1) ....  m 


In  order  to  distinguish  the  set  of  constraints  satisfied  as  equality 
constraints  at  the  solution  introduce  the  basic  set  of  constraints 


{i  | c^x*)  = 0; 


i = 1 .....  m 


Assume  that  the  gradients  of  the  constraints  in  B are  linearly 
independent  at  x . 

The  solution  of  the  nonlinear  programming  problem  has 

received  a great  deal  of  attention  in  recent  years  and  a number  of 

algorithms  have  appeared  in  the  literature.  The  accelerated  multiplier 
3 

algorithm  was  used  to  solve  the  stated  orbit  transfer  problem.  Essen- 
tially the  algorithm  consists  of  a number  of  cycles  through  the  following 
steps: 
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Step  1.  Ftnd  a point  where  7 L(x,  X)  = 0,  for  fixed  X,  by 
finding  the  unconstrained  minimum  of  an  augmented  performance 
index.  A rank  one  recursive  algorithm  is  used  for  the 
unconstrained  optimization  process. 

Step  2.  Estimate  the  Lagrange  multipliers  X,  by  minimizing 
the  error  in  the  Kuhn-Tucker  conditions.  Using  the  multipliers, 

v 

make  an  estimate  B for  the  basic  set  of  constraints  B . 

Step  l.  Extrapolate  to  find  a point  where  c^(x)  = 0,  assuming 
the  objective  function  is  quadratic  and  the  constraints  are  linear. 

| Step  4.  Reestimate  the  multipliers  X,  and  the  basis  B. 

Step  5.  Test  for  Convergence. 

Typically  two  or  three  cycles  are  required  to  obtain  the  desired  accuracy. 
For  a more  detailed  description  of  the  nonlinear  programming  algorithm, 

the  reader  is  referred  to  Ref.  3-5. 

4.  ORBIT  TRANSFER  RESULTS 

The  nonlinear  programming  algorithm  given  above  was  applied 
i to  the  stated  orbit  transfer  problem  and  Figure  2 summarizes  the  results, 

j After  coasting  in  the  park  orbit  the  first  burn  added  7820  fps  to  the 

inertial  velocity  vector  with  very  little  pitch  offset  (-1. 36°)  and  yaw 
offset  (5.  99°).  The  resulting  transfer  orbit  had  an  apogee  altitude  of 
17,  617  n mi,  while  keeping  the  perigee  altitude  and  inclination  essentially 
unchanged.  After  coasting  to  an  altitude,  h,  of  10,  838  n mi,  where  the 
argument  of  latitude,  u,  was  37.7°,  the  second  burn  was  performed.  An 
! impulsive  velocity  of  4054  fps  was  added  with  a yaw  left  of  74.  7°  and 

j a positive  pitch  angle  of  7, 21°,  Because  of  this  out  of  plane  maneuver, 

| most  of  the  required  plane  change  was  accomplished,  the  second  transfer 

orbit  having  an  inclination  of  58.  3°.  Notice  also  that  the  second  transfer 
orbit  has  an  apogee  altitude  of  23,  386  n mi  and  a perigee  altitude  of 
1,  314  n mi  and  hence  has  more  energy  than  the  desired  final  orbit. 

The  third  and  final  burn  occurred  when  the  argument  of  latitude  was 
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SECOND  TRANSFER  ORBIT 
hg=  23,386  nmi 
hp=  1,314  nmi 
i = 58.3° 

SECOND  BURN 
> V AV2  =4,054  fps 

\ A02=7-21° 

\ A^-74-7° 

\\  h =10,838  nmi 

m /u  =37.7° 


FINAL  ORBIT 

h *21,450  nmi 

tip  = 350  nmi  ^ 

i =63.4° 

u =270° 

P 


B hp  = 158.0  nmi 
JF  / i = 37.8° 

FIRST  BURN 
AV^-  7,820  fps 
^=-1.36° 
A^*5.99° 
u - 255° 


Figure  2.  Optimal  Three  Burn  Transfer,  1 = 63.4* 


121.7°  at  an  altitude  of  12,  832  n mi,  and  was  essentially  a retro  burn 
which  added  1,514  l pa  velocity  with  a yaw  attitude  of  121.2°  and  pitch 
down  of  21.  8°.  The  retro  maneuver  corrected  the  argument  of  perigee 
error  in  the  second  transfer  orbit  while  also  reducing  the  energy  of  the 
final  orbit. 

After  obtaining  results  for  the  final  orbit  with  inclination  of  63.  4°, 

a number  of  other  optimal  solutions  were  obtained  for  larger  inclinations. 

In  particular, Figure  3 presents  the  optimal  three  bum  orbit  with  an 

inclination  of  116.6°,  which  required  a plane  change  of  79.2°  from 

the  park  orbit.  All  other  constraints  in  the  original  problem  statement 

remain  unchanged.  The  first  burn  added  9686  fps  velocity  approximately 

along  the  inertial  velocity  vector  (A  0,  = -.113°,  A 4>,  - 1.080).  As  in 

the  first  case, the  resulting  transfer  orbit  remained  nearly  unchanged  in 

all  orbital  elements  except  the  apogee  altitude  which  was  increased  to 

72,  364  n mi.  After  coasting  to  an  altitude  of  65,  041  n mi,  an  additional 

3,  572  fps  was  added  to  the  velocity  in  the  direction  defined  by  A 0^  = 42.  7° 

and  = -125°.  The  resulting  transfer  orbit  had  an  inclination  of  118.4°, 

more  than  the  required  116.6°  and,  furthermore,  the  apogee  altitude  was 

7 5,000  n mi.  Note  that  for  thij  case  the  transfer  orbit  apogee  inequality 

constraint  (7)  is  satisfied  as  an  equality  constraint,  that  is^c,  (x)  is  in 

♦ ^ 
the  basic  set  of  constraints  B . The  third  burn  was  again  a retro 

maneuver,  adding  3522  fps,  in  a direction  yawed  187°  from  the  inertial 

velocity,  and  pitched  down  21. 4°. 

Figure  4 presents  a comparison  of  the  minimum  total  velocity 
increments  for  two  and  three  burn  transfers  to  a range  of  final  orbit 
inclinations.  Notice  that  there  is  a reduction  of  approximately  5.  12% 
in  the  total  AV  required  to  perform  the  mission  when  the  three  burn 
approach  is  used  for  a 26°  plane  change  (i  = 63.4°).  The  percentage 
improvement  increases  dramatically  when  the  plane  change  is  79,  2° 

(i  = 116.6°)  to  37.  5%.  It  should  be  noted  that  the  5.  12%  reduction 
in  total  AV  results  in  nearly  20%  more  payload  capability  for  typical 
IUS  configurations. 

It  is  interesting  to  observe  that  for  a plane  change  greater  than 
41°,  the  transfer  orbit  apogee  altitude  constraint  was  in  the  optimal 
basis,  whereas  for  a plane  change  less  than  41°  the  constraint  was  not  violated. 
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SECOND  BURN 
AV2  -3,572  fps 
A$2  =42.7° 


FIRST  TRANSFER  > 
ORBIT 

h = 72,364  nmi 

D 

hp=  158  nmi  — 


i = 37.5° 
w - 263.5° 


FINAL  ORBIT 

h = 21,450  nmi 
° ^ 
Hp^  350  nmi 

i *=116.6° 


SECOND  TRANSFER  ORBIT 
ha=  75,000  nmi 
\ hp=  2,692  nmi 

\ / i =118.4° 

ya)p  = 234.7° 

\ THIRD  BURN 
\ AV3=  3,522  fps 
A03=-21.4° 
A^3=  187° 


FIRST  BURN 
AV<|=  9,686  fps 
Afl^ -0.113° 
Afp  1-08° 
u = 263° 


v Optimal  Three  Burn  Transfer,  i = H*>.  6 


Figure  4.  Impulsive  Burn  Incremental  Velocities 
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The  numerical  sensitivity  of  the  sclution  in  this  region  is  apparent 
from  the  rapidly  changing  behavior  of  the  variables  in  this  region.  It 
seems  clear  that  if  the  constraints  did  not  prevent  escape,  the  optimal 
solution  would  be  to  perform  the  second  burn  at  an  infinite  altitude.  In  this 
situation  one  could  accomplish  all  of  the  plane  change  with  AV^  = 0.  This 
of  course  is  not  a realistic  solution  to  the  problem.  In  general,  for  a 
three  burn  transfer  requiring  a plane  change,  one  can  expect  that  for  a 
large  enough  plane  change  the  optimal  solution  requires  the  second 
burn  at  infinity.  The  threshold  value  for  the  specific  orbit  transfer 
considered  was  approximately  41°  plane  change. 
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EXTENSIONS  AND  GENERA LIZATIONS 


a)  Multiple  Stage  Vehicles 

The  results  presented  so  far  do  not  involve  the  vehicle 
mass  in  any  of  the  computations.  When  dealing  with  a restartable  single 
stage  vehicle, this  approach  can  be  valuable.  However,  when  a multiple 
stage  vehicle  is  considered,  the  majs  changes  that  occur  at  staging 
cannot  be  overlooked.  In  particular  one  can  consider  computing  the 
AV  for  an  N-stage  vehicle  according  to  the  equations. 


J 


A V . = g I. In 
J o j 


N 

J +£ 


<pk + V 


— F 
J + V 


JU  (Pk  + S ) - p 
k=j  * * 3 


, j = 1,  2,  . . . N, 


(17) 


where  J is  the  weight  of  the  ^tyload,  the  weight  of  the  propellant  in 
the  k th  stage,  the  weight  of  the  structure  of  the  k th  stage,  1^  the 
effective  specific  impulse  of  the  k th  stage  and  g^  the  gravitational 
constant. 

By  us;ng  the  equations  (17)  to  define  the  velocity  increments  in 
the  trajectory  simulation  described  in  Section  2a)  one  can  consider  a 
number  of  different  optimization  problems  Specifically,  one  can  pose  a 
maximum  payload  orbit  transfer  Ly  defining  the  objective  function  in  Eq. 
(1)  to  be  the  payload  J.  Sine*-  the  velocity  increments  are  defined  by 
Eq.  (17),  they  are  not  treated  as  optimisation  variables.  Instead  the 


variable  set  must  contain  the  payload  J.  When  designing  a vehicle 
one  may  also  want  the  propellant  and/or  structure  weights  and 
to  be  treated  as  variables.  Cases  have  been  run  using  each  of  the 
different  sets  of  variables,  for  preliminary  design  analysis  of  the  IUS 
vehicle.  It  should  be  clear  that  the  maximum  payload  capability 
to  a given  orbit  does  not  necessarily  result  by  flying  a minimum  £V 
orbit  transfer,  unless  the  vehicle's  weights  and  have  been 
specifically  designed  for  that  orbit. 

b)  Finite  Burn  and  Oblate  Earth  Effects 

In  the  simulation  an  impulsive  velocity  approximation 
to  a finite  burn  was  used,  and  oblate  earth  effects  were  neglected  in 
the  geopotential  model.  To  assess  the  effect  of  these  approximations, 
the  results  were  compared  with  those  obtained  from  a detailed 
simulation.  The  variables  for  the  detailed  simulation  included  pitch 
and  yaw  attitudes  at  the  beginning  of  each  burn,  co  \stant  pitch  and  yaw 
rates  during  the  burns,  the  burn  times,  and  locatior  of  each  burn.  The 
payload  difference  for  typical  missions  was  quite  small  (normally  within 
a few  percent). 

c)  Mirror  Image  Solutions 

It  should  be  observed  that  the  optimization  algorithm 
used  to  obtain  the  results  presented  only  obtains  a local  optimum  point. 
Another  local  solution  to  the  problem  doefc  exist  which  we  refer  to  as  a 
mirror  image  solution.  Because  of  the  problem  symmetry,  the  optimal 
value  of  the  objective  function  is  the  same,  although  the  location  of  the 
burns  are  different.  For  example,  in  the  mirror  image  solution  to  the 
results  presented  in  Fig.  2,  the  argument  of  tat'tude  at  ths  start  of 
the  second  burn  is  in  the  second  quadrant  rather  than  in  the  first. 

Similarly,  the  argument  of  latitude  for  the  beginning  of  the  third  burn 
is  in  the  first  quadrant  rather  than  the  second.  Although  the  objective 
function  is  the  same  between  the  mirror  image  and  the  regular  solution, 
there  may  be  some  other  criteria  for  preferring  one  to  the  other.  For 
example,  the  total  transfer  time  or  ground  track  for  one  case  may  be  more 
desirabls  than  the  othfr. 
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6. 


SUMMARY  AND  CONCLUSIONS 


The  report  outlines  a nonlinear  programming  approach  for 
obtaining  the  optimal  three  burn  solutions  to  a class  of  orbit  transfer 
problems  requiring  large  changes  in  orbital  inclination.  For  this  type 
of  problem, the  first  burn  raises  the  apogee  altitude  of  the  first  transfer 
orbit,  and  the  second  burn  is  performed  at  a high  altitude.  Since  the 
velocity  is  low,  most  of  the  plane  change  can  be  efficiently  accomplished 
by  this  burn.  The  final  burn  thus  decreases  the  energy  to  meet  the 
final  orbit  constraints.  There  can  be  a significant  AV  benefit  derived  from 
a three  burn  transfer  with  respect  to  a two  burn  transfer  as  Figure  4 
indicates.  Generalisation  of  the  approach  to  multiple  stage  vehicles  as 
as  well  as  more  detailed  trajectory  simulations  1b  discussed.  The  method 
is  currently  being  used  to  solve  a wide  variety  of  trajectory  optimization 
problems  qf  varying  complexity. 
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